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Abstract 

The high energy photoproduction of jets is being observed at the ep collider, HERA. 
It may be that the HERA centre-of-mass energy is sufficiently large that the production 
of more than one pair of jets per ep collision becomes possible, owing to the large number 
density of the probed gluons. We construct a Monte Carlo model of such multiparton 
interactions and study their effects on a wide range of physical observables. The conclusion 
is that multiple interactions could have very significant effects upon the photoproduction 
final state and that this would for example make extractions of the gluon density in the 
photon rather difficult. Total rates for the production of many (i.e. > 3) jets could provide 
direct evidence for the presence of multiple interactions, although parton showering and 
hadronization significantly affect low transverse energy jets. 



1 Introduction 



In recent years, the TEVATRON [pp) and HERA {ep) colliders have made it possible to study 
the standard theory of strong interactions (QCD) in a new regime: the regime of high parton 
densities. QCD predicts a rapid build up of very slow gluons within hadrons, which can be 
observed at high centre of mass (CM) energies, i.e. x ~ Q'^/s -C 1 where s is the CM energy 
and is the appropriate hard scale. Both the TEVATRON and HERA can study this region 
in some detail, through the production of high pt jets, heavy flavours and large-t diffractive 
scattering. At HERA, we have the additional possibility of studying this 'small-x' physics using 
deep inelastic scattering (where x = xbj). 

In this paper, we study jet production at HERA in those events where the photon is nearly 
real (photoproduction). Photoproduction of jets has been observed at HERA |l| and different 
components of the cross section identified 0. The components are easy to define in leading 
order QCD and are termed direct and resolved. In the direct sector, the photon carries all of its 
energy into the hard scatter whilst in the resolved sector, only a fraction of the photon energy 
participates in the hard subprocess. At higher orders this simple distinction between direct and 
resolved is no longer uniquely defined. However, a precise definition of the separation between 
resolved and direct processes in terms of physical observables is possible ^ allowing the effects 
of the parton distributions in the proton and photon to be disentangled [Q. This has already 
been implemented in higher order dijet cross section calculations 

The total rate is high for the photoproduction of jets and allows a study of small x phenomena 
when X ~ Ej^^ /s <^ 1 (Ej^^ is the jet transverse momentum). Naively, one might conclude that 
the 1800 GeV CM energy available at the TEVATRON would make it far superior to HERA 
(with a 7p CM energy that is typically ~ 200 GeV) regarding the study of small x physics 
through jet production. However, this is not the case. Associated with higher CM energies is 
an increase in the background of soft physics. At HERA, the asymmetric configuration of the lab 
frame (27 GeV leptons collide with 820 GeV protons) tends to boost the poorly understood 'soft' 
physics into the proton direction and hence down the 'forward' beam hole. The asymmetric 
boost of the jp system in the lab also means that jets that have been produced by small x 
partons within the proton appear in the central region in the lab frame of reference and are 
thus clearly visible in the detectors. None of these benefits are present at the TEVATRON, 
where the lab frame is also the pp CM frame. 

Due to the proliferation of low x partons (which can be inferred from the strong growth observed 
in the HERA data on the proton structure function 0), it is possible that more than one pair of 
jets can be produced per •jp collision. This multiple (parton) scattering, illustrated in figure 1, 
is expected at high enough CM energies and signals the onset of unitarization corrections to 
the simple perturbative QCD picture of 2 — >^ 2 parton scattering. Unitarization corrections are 
certainly necessary since the total cross section for inclusive jet photoproduction (calculated in 
lowest order QCD with steeply rising parton distribution functions) will ultimately exceed the 
total 7p cross sectionQ As we shall see, this apparent anomaly is resolved once we appreciate 
that the inclusive jet cross section exceeds the total cross section for jet production by a factor 
equal to the mean multiplicity of multiple interactions. We introduce an eikonal model which 

^Multiple scattering is only one manifestation of the unitarization corrections that are expected to occur at 
high energies. The steepness of the parton densities will eventually be tamed by additional corrections. We 
ignore such effects in this paper. 
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Figure 1: An example of a multiple scattering in a 7p collision. 

assumes that individual hard scatters are uncorrelated. This allows us to model the rate of 
multiple interactions and study their effects upon the hadronic final state. For the latter study, 
it is most convenient to subject our model to a Monte Carlo simulation and this will allow us 
to make full use of the available and forthcoming HERA data. 

In section 0, we describe the eikonal model of refs. 0, ||, Q and show how it predicts a significant 
rate for the production of multijet events at HERA energies. Also described here is the imple- 
mentation of the formalism within the HERWIG Monte Carlo package By integrating our 



formalism within HERWIG, we can make realistic studies of the final state which include the 
effects of parton showering and hadronization. In section § a number of key photoproduction 
distributions are presented for our default model, and compared to the results obtained with- 
out multiple interactions. In sections 4 and 5, we examine the effects of using different parton 
distribution functions and of variations on the default model. In section 6, comparisons are 
made to available HERA data. We show that multiple scattering can be expected to make a 
clean extraction of parton distribution functions in the photon rather difficult. 



2 Why Multiple Interactions at HERA? 

We are interested in jet production in 7p reactions and at the level of the hard subprocess we 
assume that this can be approximated by the lowest order matrix element for 2 — >■ 2 parton 
scattering with final state partons produced with a transverse momentum, px > Pt™ 3> Aqcd- 
Jets can be produced directly via the 7-parton hard subprocess or indirectly via partons from the 
resolved photon scattering with partons in the proton. These resolved partons can be generated 
either non-perturbatively (i.e. the 7 splits into a large size qq pair) or perturbatively (i.e. via 
perturbative evolution of a small size qq pair). Let us start by considering a interaction at 
some fixed centre-of-mass energy, s^p. In the CM frame we think of the proton and resolved 
photon as Lorentz contracted 'parton pancakes' colliding at some impact parameter, b. The 
mean number of jet pairs produced in this resolved-7-p interaction is then 

{n{b, S)) = ^partons ® (l) 

where /^partons is the parton luminosity and cth is the cross section for a pair of partons to 
produce a pair of jets (i.e. partons with px > p™°). The direct photon interaction generates 
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only a single pair of partons and we ignore it for now (it will be included as a separate hard 
subprocess in the Monte Carlo simulation). 

The convolution is because the parton cross section depends upon the parton energies. More 
specifically, 

'^'^partons A{l))Ti^{^X^^Tlp{Xj^dx.ydXp (2) 

where ni{x,i) is the number density of partons in hadron i which carry a fraction x.i of the hadron 
energy. For ease of notation we do not distinguish between parton types and have ignored any 
scale dependence of the number densities. A{h) is a function which specifies the distribution of 
partons in impact parameter. It must satisfy 

Tldh'^A{h) = 1 

in order that the parton luminosity integrated over all space is simply the product of the parton 
number densities. Factorizing the b dependence like this is an assumption. In particular we 
do not contemplate QCD effects which would spoil this, e.g. perhaps leading to 'hot spots' of 
partons. For the proton, the number density is none other than the proton parton density, i.e. 
np{xp) = fp{xp). The number density of partons given that the photon is resolved is related 
to the photon parton density by a factor of ~ ttem, i-e. n^{x^) ~ f^{x^)/aem- In our model, we 
ignore multiple scattering which arises due to interactions of the small size fluctuations of the 
incoming photon and therefore assume that the photon interacts just like a hadron. The large 
size hadronic fluctuations of the photon can be modelled assuming vector meson dominance 
and allow us to estimate that n-y(x-y) = /y(x-y)/Pres where Pres = KATraem/ fp- The sum over 
hadronic fluctuations determines k. Unless otherwise stated, we assume p-dominance, i.e. k — 1 
and Pres ~ 1/300. 

Thus, after performing the convolution, we can write: 

{n{b,s,p)) = ^a^^{s,p), (3) 

-'res 

where cr^^ls^p) is the inclusive cross section for 7p — > jets. Restoring the parton indices, it is 
given by 

In order to investigate further the structure of events containing multiple interactions we need 
to know the probability distribution for having m (and only m) scatters in a given resolved-7-p 
event, P^. In order to do this we assume that the separate scatters are uncorrelated, i.e. they 
obey Poissonian statistics. Thus 

Pn. = exp(-(n(6, (5) 

This formula is central to the Monte Carlo implementation in HERWIG. 
We can now ask for the total cross section for 7p — > partons with px > 



(^h{s^p) = TrPres j db'^ ^ Pm 

= TrPres y c?6^[l - exp(-(n(6,s^p)))]. (6) 



3 



Since the total inclusive cross section (o")^'^) counts all jet pairs (even ones which occur in the 
same event) we expect it to be larger than au by a factor equal to the mean number of multiple 
interactions per event (i.e. averaged over impact parameter). This is easy to see. Let (n(s^p)) 
be the average number of jet pairs produced in resolved-7-p events which contain at least one 
pair of jets, then 



(n(s^p)) 



Jdb^{n{b, s^p)) 
Jdb^[l -exp(-(n(6, s^p)))] 

Note that an must always be less than the total 7p cross section, whereas a^^^ need not be. 



(7) 



Jet cross sections at HERA are sensitive to the proton parton density down to small momentum 
fractions, x ~ 10~^. As shown by the measurement of the structure function in deep inelastic 
scattering at HERA, this is the region of high parton density (which rises as at ~ 

10 GeV^ 0). For pf"" ~ 1 - 2 GeV, one expects that aj^"(l TeV) ~ 1 mb, which is much 
larger than the total jp cross section, which is ~ 200 fih according to the successful Regge 

This suggests the possibility that {n{s^p)) is 



pole model of Donnachie and Landshoff p2| 
significantly greater than 1. 

To complete the description of our model, we must now specify how the partons are distributed 
in impact parameter space. As a first approximation, we ignore pQCD effects and assume that 
A{b) can be written as a convolution of form factor distributions, i.e. starting from 



Gp{h) 



d'^k exp(k-b) 
(27r)2 (1 + kV/i2)2 

d'^k exp(k ■ b) 



(27r)2 1 + kVz/2 
with = 0.71 GeV^ and u'^ = 0.47 GeV^ we can write 



A{b) = J d^h' G'p(b')G'^(b - b') 



(8) 



(9) 



The integral then yields 



271 /i2 



/i 



:{Ko{ub) ~ Koifib)) - ^K.ifib) 



(10) 



The formalism just described has been discussed extensively in the literature (see ref. |T3|, |T 



but largely in the context of total cross sections. The expected rates for multi-jet processes 
have also been studied [|l^. Owing to the complexity of the final state, in which many jet 



pairs may be produced, previous studies have been confined to rather inclusive properties of 
the model. To really test for the existence of multiple interactions however, it is necessary to 
confront the vast amount of data that HERA is making available with more detailed properties 
of the final state and to this end we have constructed a Monte Carlo simulation program. In 
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ref.|16] some results of an earlier version of this simulation were presented. In this paper we 
present a more comprehensive study of the effects of this model on the photoproduction final 
state. The program is written as an extension to the HERWIG [|I0| package and can also be 



used to simulate the effects of multiple interaction in other processes, such as 77 interactions in 



e'^e collisions |]T^. There is a degree of overlap with the multiple interaction formalism that 



is included in the PYTHIA program 



The effect of eikonalization on the total 7-p cross section is shown in figure |^. In figure 03 
this is translated into the effect on the differential ep cross section dcjep/ dy where s^p = AyE^Ep 
{Ee is the incoming electron energy and Ep the incoming proton energy). The p™™ is taken 
to be 3 GeV and Pros = 1/300. The parton distribution sets used are the GS2 set [jl9| for the 
photon and the GRV [|20| set for the proton. Clearly (for this choice of parameters) the effect of 
multiple scattering is to significantly lower the cross section, with significance increasing with 
increasing s^p (and therefore y). This is reflected by a significant probability for more than 
one hard interaction in a 7p event, as shown in figure |^c, where it can be seen that, depending 
upon y, as many as 10-15% of 7p events contain more than one parton-parton scattering. This 
is expressed as a differential cross section for N and only scatters in figure |^d. Note that 
the hard cross section for ^/s-yp ~ 200 GeV represents around 20% of the total 7-p cross section 
flTTl. 



2.1 Monte Carlo Implementation 

The cross sections illustrated in figure ^ indicate that the effects of multiple interactions might 
be observable at HERA. However, the experiments do not measure jet cross sections as low as 
3 GeV. The lowest i^^?* cut used on measured cross sections so far has been 6 GeV. The 
smearing of jet energies due to parton showering and hadronization effects mean that a p™™ 
cut of around 2.5 GeV is often used by the experiments in Monte Carlo simulations to obtain a 
good description of the data at higher Ej^^ cuts. The pile-up of more than one hard scatter in an 
event can be expected to increase these migrations, and in order to study the effect of multiple 
interactions on real jets rather than partons, it is necessary to include realistic hadronization 
and parton showering models in conjunction with the eikonal model described above. To this 
end, the HERWIG Monte Carlo program was adapted to allow the generation of multiple hard 
subprocesses in a single hadron-hadron, photon-hadron or photon-photon event. 

In the default version of HERWIG, event generation begins with the hard process, QCD dijet 
production in our case. This is generated according to the lowest order cross section using 
whichever parton distributions have been selected by the user. Associated with the hard process 
is a hard scale, Q (~ E^:^ for dijet production), and the cross section for additional partons to 
be radiated in the same event is enhanced by powers of logQ/Qo at each order in as, where 
Qo is a lower cutoff defining the resolvability of partons. In order to fully describe the final 
state, Qo is required to be around the typical hadronic scale, Qo ~ 1 GeV, and the logarithmic 
terms are large. Therefore such terms must be summed to all orders in a^. This summation 
can be implemented as a probabilistic evolution of the external (both incoming and outgoing) 
partons. Any parton shower algorithm that implements the DGLAP splitting functions 
correctly resums the single logarithmic terms associated with the collinear emission that drives 
fragmentation function and parton distribution function evolution, but to fully model the final 
state, single and double logarithmic terms associated with soft emission must also be summed 
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Figure 2: a) an as a function of s^p, the squared cm. energy of the photon-proton system, 
b) da'j^/dy as a function of y. In a) and b) the dashed hnes show the resuh with no muhiple 
interactions, the sohd hnes show the resuh of including multiple interactions and the points 
indicate the cross sections actually generated by the Monte Carlo program (see text). In c) the 
probability of N and only N scatters as a function of y is shown. In d) the cross section for 
and only N scatters as a function of y is shown. 
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to all orders. This is done by coherent parton shower algorithms , which also implement the 



DGLAP splitting functions, but with specific choices for the definitions of evolution scale and 
parton resolution criterion. For final-state partons, the algorithm evolves a single hard parton to 
many softer ones. For initial-state partons, the evolution is 'backwards', from the hard process 
to the incoming hadron or photon, evolving a single incoming parton to another at larger x, 
together with many additional outgoing partons. The backward evolution algorithm [2^ ensures 



that at each stage of the evolution, the distribution of parton x values agrees with that of the 
input distribution function. For the photon, the inhomogeneous term in the evolution equation 
appears as an additional vertex in the backward evolution, giving the possibility that it will 
terminate 'anomalously', with the incoming photon giving all of its energy to the hard process 
and its associated parton shower. 

Finally, after perturbative parton evolution, a non-perturbative model of the transition from 
partons to hadrons is invoked. In HERWIG, a simple model that makes use of the preconfine- 
ment property of QCD is used. At leading order in the number of colours, the structure of a 
parton shower is such that each parton ends up close, in both momentum-space and real-space, 
to a parton that carries the opposite colour quantum numbers. The model therefore combines 
each such colour-connected pair of partons to form a single colourless cluster. This is assumed 
to decay directly to hadrons, with the choice of hadron species controlled purely by phase-space 
arguments This model therefore has few adjustable parameters, which are already rather 
well constrained by data from e"'"e~ annihilation. For processes with incoming partons, there is 
an additional freedom, of how to model the break-up of the incoming hadron or photon after 
one parton is taken from it. In HERWIG, this is done by requiring that the backward evolu- 
tion results in a valence parton (u, d and s quarks and antiquarks are all considered valence 
partons of a photon). If the perturbative evolution did not end anomalously, then additional 
emission is generated below Qo, but still according to the perturbative distributions, to obtain 
a valence parton. Once the 4-momentum of this valence parton is known, energy-momentum 
conservation is used to calculate the momentum of the remaining hadron or photon remnant. 
For photons, the fiavour of the remnant is simply that of the antiparticle of the valence par- 
ton, while for protons, the remaining diquark is taken to be a single anticoloured parton. The 
cluster containing the remnant parton is then hadronized like any other. Clearly this step is 
unnecessary in events where the backward evolution terminated anomalously, since there is no 
remnant parton in that case. 

When incorporating our multiple interactions model into the Monte Carlo event generator 
several additional problems arise from the fact that the model assumes that different scatterings 
are uncorrelated, while in the event generator this is not possible. The first is that of energy 
conservation. In the analytical model, there is a non-zero (though small) cross section for 
events in which the total energy participating in hard scatters is greater than the energy of the 
incoming particles. This is exacerbated by the backward evolution, since additional energy is 
radiated away in the additional partons. We impose global energy-momentum conservation by 
adding the simplest possible correlation between scatters: a 9-function in the cross section to 
produce m and only m scatters requiring their total energy to be less than the available energy. 
This will reduce the amount of eikonalization of the cross section relative to the analytical 
model. When generating events, this is implemented by calculating the total cross section and 
probabilities for exactly m scatters using the analytical model, and generating hard scatters 
according to these probabilities. Scatters in which energy-momentum cannot be conserved are 
then vetoed (although the event itself is kept, with all the other scatters). At the end of event 
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generation, the number of vetoed scatters is used to make a revised estimate of the total cross 
section. 

In the analytical model, only events in which the photon is resolved are eikonalized, while 
direct photon events are not. However, events in which the backward evolution terminated 
anomalously should also be called direct in this context (they correspond to the small size 
fluctuations of the photon) and should not be eikonalized. Therefore, if the backward evolution 
does terminate anomalously, no further scatters are allowed in that event. This is treated just 
like the vetoing from energy conservation so also revises the final estimate of the total cross 
section. Since the separation into the anomalous and resolved parts of the distribution function 
is model dependent, we provide an option that enables multiple interactions in anomalous 
events, so that one can gauge the relative importance of this feature of the model. 

These revised cross sections are shown as data points on figure |^a and b. 

The simple model used in HERWIG for the hadron and photon remnants is no longer appropri- 
ate when there is more than one scatter. This is because once the incoming colourless hadron 
has been replaced by an outgoing coloured remnant by the first interaction, the same procedure 
cannot be iterated for subsequent interactions. Instead, we label the remnant as a new kind 
of incoming hadron for the other interactions, which has identical properties to the original 
hadron, except that gluons are labelled as its valence partons. Thus, if the backward evolution 
does not result in a gluon, an additional emission is forced to produce one. The outgoing 
coloured remnant from this interaction is the same as that from the first, but with reduced 
momentum. This procedure can then be iterated as many times as required for all subsequent 
interactions. By thus modifying only the part of the model that connects the hard processes to 
the incoming hadron, and not the backward evolution itself, we ensure the predictivity of the 
model, using the parameters fixed to other reactions. 



3 Effect on the Hadronic Final State 



The jet properties of the final state are of great importance in understanding the underlying 
QCD processes, and the possible influence of multiple scattering needs to be carefully examined. 
For example, jet rates in photoproduction may well lead to important information regarding 
the gluon content of the photon and the proton [|, In this section we examine the general 



effects of multiple scattering on the properties of the final state measurable at HERA. In order 
to do this, we take a default set of parameters chosen where appropriate to reflect published 
data from the HERA experiments and the best available theoretical estimates. These choices 
are as follows: 



-Pres = 1/300: This is motivated by assuming p-dominance. 

p™™ = 3.0 GeV: This is consistent with the values typically used by the experiments to 
describe their data at Ej^^ > 6 GeV. 

The GRV parton distribution set is used for the proton [^. This set is the result of a 
global analysis including the HERA measurements of the proton structure function F2 H. 
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• The GS2 parton distribution set is used for the photon . The gluon distribution in the 
photon is poorly constrained over the region studied (the desire to constrain the gluon 
is a major motivation for measuring high-i^^T photoproduction cross sections at HERA) 
and there is little reason to favour any one out of several available parameterizations. 

• A cut on the 7p CM energy was made, 114 GeV < ^Js^p < 295 GeV, similar to those 
usually made by the experiments. 

These choices define what we call our 'default' model. In the subsequent sections the sensitivity 
to various choices made here will be examined and discussed. Events were generated both with 
and without multiple interactions. All other parameters of the Monte Carlo model were left 
at their default values, which are tuned to e"'"e~ data and give a good description of the direct 
component of photoproduction at HERA 0. 

The initial expectation is that each pair of jets will be produced back-to-back azimuthally with 
balancing transverse momenta. For two-to-two parton scattering in leading order QCD energy 
and momentum conservation gives 



LO X]partons(-^ ~l~ P2)parton LO X]partons(-^ Pz)psxton 



LO ^partonsV" ' ^'zyparion ^partonsV" t^zjpa.rmii ('\-\\ 



The proton direction defines the positive z-axis, E^. and are the initial lepton and proton 
energies and the sum is over the two final state partons. For direct photon events, xh^ = 1. 
Since it is not possible to measure partons an observable has been defined by ZEUS in terms 
of jets which is analogous to xl^'^. This observable, called x"^^, is the fraction of the photon's 
momentum which emerges in the two highest E"^* jets. The explicit definition is. 



~ 2yE, ' ^^^^ 



where now the sum runs over the two jets of highest i?^^*. In the x"*^*^ distribution thus ob- 
tained, the LO direct and resolved processes populate different regions, with the direct processes 
concentrated at high values of x'^'^. The peak arising from the direct contribution will not nec- 
essarily lie exactly at x'^^ = 1 due to higher order effects and/or hadronization, but will still 
correspond to the kinematic region where most or all of the energy of the photon participated 
in the hardest scatter. Making a cut on x"^^ provides a workable definition of direct and re- 
solved events. Recall that beyond lowest order, such a separation is ambiguous. We introduce 
a similar definition for the proton, i.e. 



obs 



2E„ 



(13) 



In order to make our study as realistic as possible, we have performed jet finding on the hadronic 
final state using a cone algorithm. Unless stated explicitly otherwise, the cone radius used is 
R = 1 {as has been used by both HERA experiments so far), and jets have Et > 6 GeV 
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and pseudorapidity 77 = — ln(tan^/2) < 2. These cuts place our jets well within the region 
observable by the HERA detectors. 

Multiple parton scattering is expected to affect jet rates in two ways. Firstly, and most obvi- 
ously, the average number of jets per event should be increased when partons from secondary 
hard scatters are of sufficiently high pj- to give jets in their own right. In addition, there is 
also the possibility that extra hard scatters can influence the observation of jets even when no 
parton from the secondary scatters is of high enough to produce an observable jet. Lower 
Pt secondary scatters produce extra transverse energy in the event which can contribute to 
the pedestal energy underneath other jets in the event. In this way, the inclusion of multiple 
scattering can be thought of as an attempt to extend perturbative QCD in order to calculate 
some part of the so called 'soft underlying event', which is included in many simulations as 
a parameterized extrapolation of existing data. The first of these effects will only increase 
multi-jet cross sections, whilst leaving the inclusive jet cross section unchanged. The second 
effect will increase both inclusive and multi-jet cross sections. In figure |, where inclusive jet 
(figure ^,b) and dijet (figure |^c,d) cross sections are shown (as a function of jet rapidity and 
mean jet rapidity respectively), it can be seen that the effect of multiple interactions is indeed 
more significant for the dijet cross sections than the inclusive jet cross sections. In all cases the 
effect is more significant at lower E^r^^. The direct photon contribution generated by HERWIG 
is included. 

In table |l| the calculated rates for one-, two-, three- and four-jet events in the 'default model' 
are presented. The cross section for m jets is defined here as the production cross section 
of events containing m and only m jets. There is an enhancement in all jet rates above the 
single jet rate and the significance rises with increasing jet multiplicity. Single jet events arise 
primarily from configurations where the second jet is forward of the rj cut. Without considering 
the extra transverse energy generated by multiple scattering, we would expect the sum of the 
jet cross sections obtained including multiple scattering to be lower than the corresponding 
sum obtained without multiple scattering (and the latter should be concentrated more at lower 
jet rates), by a factor equal to the mean multiplicity of jet pairs. This effect is mostly washed 
out by the large effect of the 'underlying event' generated by the secondary scatters. In this 
table only the contribution from LO resolved photon processes is included. Since the parton 
shower model without multiple scattering only generates multi-jet cross sections as a leading 
logarithmic dressing of the dijet cross section, it is only accurate in the region in which two jets 
are significantly harder than all the others. Therefore the predicted rates should only be taken 
as a rough indication. Nevertheless, since the logarithmically enhanced terms are included, the 
neglected terms are genuinely suppressed by powers of a^, without additional enhancement, 
whereas the eikonal expansion is an expansion (summed to all orders) in the large parameter 
agN"^ (where A^^ is the product of the number densities). We therefore expect the enhancement 
of higher jet rates to remain significant. 

These results imply that, even with these rather strict cuts on the jets, around 500 four- 
jet events (where the jets are easily observable in the detector) generated from multiple hard 
interactions should be present in the data already taken by each of the HERA detectors (around 
10 pb-i by the end of 1995). 

Figure ^ shows the x'^'^ and x'^'^ distributions, again including the direct contribution generated 
using HERWIG (which can be seen peaking at high values of x'^'^). The inclusion of multiple 
scattering raises the cross section in the low x'^^ region without affecting the high x^^ region, 
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Figure 3: Inclusive jet cross sections: a) da/drj and b) da/dEr- Dijet cross sections: c) da/drj 
and d) da/dEx- In all cases, i?^* > 6 GeV. The solid lines show the distributions when multiple 
scattering is included, the dashed lines show the distributions when no multiple scattering is 
allowed. In the E^r^^ plots, the statistical errors are indicated on the distribution for which 
multiple interactions are included. The errors on the 'no multiple interactions' distribution are 
similar. 
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Cross section (nb) 




No multiple scatters 


With multiple scatters 


One jet 


90.8 ±0.6 


89.9 ±0.8 


Two jet 


17.0 ±0.3 


19.5 ±0.4 


Three jet 


0.73 ±0.06 


1.15 ±0.09 


Four jet 


0.02 ±0.01 


0.05 ±0.02 



Table 1: Resolved photon cross sections for multijet events. Jets have E^: > 6 GeV, rj < 2, 
R = 1 and the proton and photon parton distribution sets used are GRV and GS2 respectively. 
The errors shown are statistical. 
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Figure 4: a) da/dx"^^ b) da/dxf"^ The sohd histo grams show the distributions when multiple 
scattering is included, the broken histograms show the distributions when multiple scattering 
is neglected. 
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Figure 5: Transverse energy flow as a function of Arj integrated over | A0| < 1.0 for jets with 
Et> Q GeV in the range -2.0 < r/ < 2.0 and x°^^ > 0.75 (a) and x°^' < 0.75 (b). The sohd 
hues show the distributions when multiple scattering is included, the dashed lines show the 
distributions when no multiple scattering is allowed. 

which is dominated by direct photoproduction. In the (uncorrected) x"^^ distributions presented 
in refs. p, ^ there is an excess of data over the standard simulations in just this region. Multiple 
interactions may have a role to play in resolving this discrepancy. The effect of multiple 
interactions can be seen to be largest for the Xp values in the middle of the available range. 
This is easy to understand since the lower Xp values are strongly correlated to large x^ values 
(where direct photon interactions are dominant). 

Figure ^ shows how multiple scattering affects the profiles of jets in the ri-(j) plane. In particular 
we look at jets in both the high and low x'^'^ regions. 

The effect of this model for multiple interactions on the jet profiles is significant but fairly 
small. In the high x'^'^ region the effect is negligible. In the low x*^^ case, there is a slight 
broadening of the jet and an increase in the pedestal of energy flow around the jet by around 
300 MeV per unit of rj on the forward (proton) side of the jet and by around 100 MeV on the 
rear (photon) side. In measured energy flows at HERA, an excess of energy flow in the forward 
region has been observed and shown to occur principally at low values of x*^^ [Q. None 
of the jet profiles presented by the HERA experiments are corrected for detector effects, and 
so smearing of jet and particle energies makes direct comparison impossible at present. One 
might take the view that the jet profile data could be used to tune p™™. This is not a very 
sensible thing to do, since (as we shall see) uncertainties in the photon gluon density can easily 
be traded off against the value of 
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Cross section (nb) 




No multiple scatters 


With multiple scatters 


One jet 


93.2 ± 1.0 


124.4 ± 1.2 


Two jet 


15.0 ±0.4 


28.3 ±0.6 


Three jet 


0.62 ±0.08 


2.13 ±0.16 


Four jet 


0.01 ±0.01 


0.15 ±0.04 



Table 2: As Table 1, but using LACl parton distribution functions for the photon instead of 
GS2. 

4 Sensitivity to Parton Distributions 

The fact that jet cross sections, and in particular dijet cross sections, are affected by the 
presence of multiple hard interactions means that these cross section may be misleading if used 
to distinguish between different parton distributions in the proton and/or photon without giving 
due consideration to the effects of multiple interactions. In figure § we show the x*^^ distribution 
and the low-x"^*^ jet profiles obtained using the GS2, LACl and GRV parton distribution sets 
for the photon. 

In general, for this kinematic region, the effect of multiple scattering is greater than or similar 
to the changes in the cross sections produced by using different parton distributions. This 
is particularly notable in the jet profiles, which in the absence of multiple interactions are 
insensitive to the photon parton distribution, but which depend upon it significantly when 
multiple interactions are allowed. The effect of multiple interactions is most pronounced if the 
LACl parton distribution set is used. This is due to the high density of low-x^ gluons present 
in the parameterization. In table ^ we show the rates of multijet events for the LACl photon 
parton distribution set, where it can be seen that the multijet rates are massively enhanced by 
multiple interactions. 

5 Variations on the Default Model 

Despite the fact that it is a reasonably simple and well defined procedure, there are of course 
several uncertainties and free parameters in the model. These include the values of p™", Pj-^s, 
and u^. The parameter which has the strongest effect of the cross sections is the value of p™'^. 
Lowering p™™ to 2 GeV dramatically increases the effect of multiple interactions, as shown by 
the effect on the jet profiles and and a;"*^*^ distribution in figures |^a and c. The energy flow 
around the jet axis and the cross section at low x?l^^ increase when p™™ is lowered. This occurs 
whether or not multiple interactions are included, but in the presence of multiple interactions 
the effect is much greater. 

The value of P^es used in our 'default model' is around the lowest sensible value. Increasing it 
from 1/300 to 1/150 decreases the effect of multiple scattering since we now correspondingly 
reduce the number density of resolved partons. Since Pres is the probability that a photon is 
resolved, its value is closely related to the photon's parton distribution functions. Indeed, if one 
fixed the parton densities in a resolved photon and only varied Pres, there would be no variation 
in the effect of multiple interactions. However, the biggest differences between current sets is 
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Figure 6: Effect of different photon parton distribution sets: a) and b) show transverse energy 
flow as a function of Arj integrated over |A0| < 1.0 for jets with Et > 6 GeV in the range 
—2.0 < Tj < 2.0 and x'^^ < 0.75. The solid and dashed lines show the distributions obtained 
using the GS2 photon parton densities with and without multiple interaction respectively. In 
a) the dotted and dash-dotted lines show the distributions obtained using the GRV photon 
parton densities with and without multiple interactions respectively, and in b) they show the 
same distributions obtained using the LACl photon parton densities. Figures c) and d) show 
the a;°^^ cross section. The lines in c) have the same meaning as in a), and those in d) have the 
same meaning as in b). In all cases the proton parton distribution set is taken from GRV. 
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in the distribution functions within a resolved photon, rather than in the value of Pres used. 
Therefore most of the differences shown in section ^ would remain even if P^es was adjusted to 
the values in the distribution functions. 

Changing the value of from 0.47 GeV^ to 1.5 GeV^ makes the photon more compact, i.e. it 
increases the resolved-7 number densities at small impact parameters whilst decreasing them at 
large impact parameters. Since peripheral 7p collisions are rarer, the net result is to increases 
the effect of multiple interactions for each resolved-7-p interaction but to make the cross section 
for such reactions smaller. The effects of these changes on the jet profiles and x°^^ distribution 
are shown in figures |^ and d. 



6 Comparison with Data 



Both ZEUS and HI ||2^ have published inclusive jet cross sections for photoproduction at 



HERA, and ZEUS has also published dijet cross sections 0. Our default multiple interactions 
model is compared to a selection of this available data in figure |^. 

The dijet cross sections are measured as a function of the average ri of the two jets if]) for 
direct and resolved photoproduction (defined by a cut at x'^^^ = 0.75) for Ej^^ > 6 GeV, 
P^ < 4 GeV^,0.2 < y < 0.8 and \Ar]\ < 0.5 (P^ is the photon virtuality). There is a 
correlated uncertainty due to the calorimeter energy scale of around 20% in the direct and 
25% in the resolved data, which here we have added in quadrature to the other systematic 
errors. The agreement with the direct dijet cross section is good in all cases, confirming that 
the separation based upon x?^^ removes to a large extent any sensitivity of this cross section 
to the photon structure (and hence multiple interactions). In the case of the resolved dijet 
cross section, the calculations are in general too low, as was also the case for the analytic 
LO QCD calculations presented in [Q. Multiple interactions raise the calculated cross section 
towards the data and, given the large systematic errors on the data (in particular the overall 
normalization uncertainty) perhaps the disagreement is not so significant yet. A major source 
of the systematic errors in the data is the discrepancy in the jet profiles between data and 
Monte Carlo models, it is to be hoped that improvements in the simulations (of which the 
inclusion of multiple interactions is an example) will contribute to the reduction of these errors 
in future measurements. 

The inclusive single jet measurements of HI and ZEUS are made over different y ranges (0.2 < 
y < 0.85 for ZEUS and 0.25 < y < 0.7 for HI), different P^ ranges (P^ < 4 GeV^ for ZEUS 
and P2 < 0.01 for HI) and different E^^' ranges (E^^* > 8GeV for ZEUS and E^^ > 7GeV for 
HI). There is a correlated uncertainty due to the calorimeter energy scale of around 20% in 
the ZEUS inclusive jet data, and a 25% overall normalization uncertainty in the HI data. In 
both cases we have here added these in quadrature to the other systematic errors. In addition 
it is worth noting here that since our model uses only the LO matrix elements plus parton 
showers, there is an expected uncertainty in the calculations due principally to the choice of 
scale. In the calculations available for these processes, the NLO corrections have been as high 
as 40% |2^ although, as mentioned earher, the parton shower would be expected to reproduce 



at least part of this correction. Since the ZEUS data are collected up to larger values of y, 
multiple interactions are expected to have a larger effect and this is seen by comparing the solid 
and dotted curves in figure 8c with the corresponding curves in figure 8d. The broken lines, 
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Figure 7: Sensitivity to model parameters: a) and b) show transverse energy flow as a function 
of Ar; integrated over |A0| < 1.0 for jets with Et > 6 GeV in the range —2.0 < 77 < 2.0 and 
^obs < 0.75. In all cases, the sohd and dashed lines show the distributions obtained with and 
without multiple interactions respectively in the default model). In a) the dotted and dash- 
dotted lines show the distributions obtained with and without multiple interactions respectively, 
when p™" is lowered from 3 GeV to 2 GeV. In b) the dotted line shows the distribution obtained 
with multiple interactions when the photon radius i/^ is changed to 1.5 GeV^ and the dash- 
dotted lines show that obtained with Pres — 1/150. Figures c) and d) show the cross 
section. The lines in c) have the same meaning as in a), and those in d) have the same meaning 
as in b). 
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Figure 8: ZEUS data compared with HERWIG. a) and b) show the differential dijet cross 
sections for direct and resolved photon events respectively (where the separation is defined by 
a cut on x'^^) by the ZEUS collaboration [Q. c) and d) shows the differential inclusive jet 
cross section as a function of 77-'°* as measured by the ZEUS |^5[ and HI collaborations 
respectively. The inner error bars indicate the statistical uncertainty, and the outer one the 
statistical and systematic errors added in quadrature. In all cases the solid histogram shows 
the calculation of the default multiple interaction model, the dotted line shows the calculation 
using the default but using LACl instead of GS2 for the parton distribution set, and the dashed 
line shows the result of HERWIG with no multiple interactions. 
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corresponding to the HERWIG predictions without multiple scattering are very similar (the 
effect of increasing the range in the ZEUS data is compensated by the lower E}^^ cut of HI). 
Given the uncertainties mentioned above, we should be cautious in drawing any conclusions 
from these data. However, there is a suggestion - at the (1 ~ 2)a level - that the HERA data 
may be incompatible, in that HI inclusive jets agree better with LACl distribution (or very 
significant multiple interaction effects) whilst the ZEUS inclusive jet data lie closer to the result 
of our default multiple interaction model. 

Bearing in mind the above discussions, there are signs that including multiple interactions can 
bring the model into better agreement with existing data in those regions where discrepancies 
exist (i.e. the forward region in both dijet and inclusive jet production), whilst having a small 
effect on those regions (i.e. direct dijet, low rj^'^^ inclusive) where the agreement is already good. 



7 Conclusions 

We have performed a detailed study of the unitarization corrections which are expected to 
appear in high energy interactions. We used an eikonal approach to model their effects 
which, in the perturbative domain, manifest themselves through the appearance of multiple 
parton scattering. By performing a Monte Carlo simulation of our model we have been able to 
make a study of the detailed properties of the hadronic final state which enables us to make 
use of the vast amount of HERA data on photoproduction. We find that, for partons produced 
with pt > 3 GeV, over 4% of events can be expected to contain more than one hard process at 
typical HERA energies. 

For reasonable experimental cuts the effect of multiple parton interactions on measured jet cross 
sections could be significant. The size of the effect is expected to depend strongly upon the 
parton distributions in the photon and can be as high as 100% in some regions for dijet cross 
sections already measured at HERA. The effect is even more significant for higher jet rates, 
leading to an overall enhancement of a factor of up to around five in the four jet rate (with all 
four jets contained in a detector). These events should provide a means of unambiguously dis- 
covering or ruling out some multiple interaction models in the near future, although untangling 
a clear signature in a realistic multi-hadron final state remains something of a challenge. 

The model considered here also leads to the generation of a 'semi-hard' underlying event. This 
affects the energy flow forward of the jet direction, as do other models which are already being 



used to to describe HERA data [|I8], where the energy flow in the forward (proton) direction 
is poorly reproduced in simulations which do not include multiple interactions. 

Collectively, the effects of multiple interactions may well make the extraction of photon parton 
distribution functions (especially the gluon) rather difficult. 

The Monte Carlo program used here runs in conjunction with HERWIG version 5.8 and is 
available from the authors. Further information can be obtained from the world wide web page 
http : / /suryal 1 . cern . ch/ users/ seymour/herwig/ 
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